AN EXPERIMENTAL INVESTIGATION 
OF THE EFFECT OF ADVERSE PRESSURE GRADIENT 
ON VORTEX BREAKDOWN 


Gale Edward Treiber 





AN EXPERIMENTAL INVESTIGATION 
| OF THE EFFECT OF ADVERSE PRESSURE GRADIENT 
* ON VORTEX BREAKDOWN 


by 


Gale Edward Treiber 


|Thesis Advisor: 
ee aaa GEIST EER 


Approved for public release; distribution unlimited. 








at Sarpkaya | 


1135104 





an 4 a 


Pa 
OAreh a. wr a ea 


An Experimental Investigation 
of the Effect of Adverse Pressure Gradient 
on Vortex Breakdown 


by 
Cale Edward) iiimned ber 
Lieutenant Commander, United States Navy 


B.S. , United States Naval Academy, 1964 
M.S., United States Naval Postgraduate School, 1971 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1973 





L Kee 
Na aval , Postgraduate School 
yvic Yonterey, California 93940 


ABSTRACT 


The results of an experimental investigation of the effect 
of adverse pressure gradient on the vortex breakdown phenom- 
enon in a diverging tube are presented. Adverse pressure 
~gradient was found to be as significant in determining the 
breakdown position as were the previously known parameters, 
namely, the Reynolds and circulation numbers. It was found 
that an increase in any of the three parameters serves to 
move the breakdown upstream, toward the origin of the vortex. 
It was further found that beyond some range, an increase in 
adverse pressure gradient causes boundary-layer separation 
in the diverging tube. Separation restricts the swirling 
flow and thus limits the effective adverse pressure gradient 
aACwiIniP On mune VOGUE x.) Davasare presenved TO aiiiusitrave 
these effects. An approximate momentum analysis for pre- 
Gicting the vortex breakdown position was carried out. The 
results were found to be only marginally satisfactory. In 
spite of that, however, the method appears to hold much 


promise as a breakdown predictor. 
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I. INTRODUCTION 


Vortex breakdown is an abrupt change in the etrueture 
of the core of a swirling flow. It was first noted by Peck- 
ham and Atkinson in 1957 during their studies of air flow 
over delta wings at large angles of attack [1]. Since then 
vortex breakdowns have also been observed in straight or 
diverging cylindrical tubes, cyclone separators, behind tor- 
pedos, and in both nozzels and diffusers. Vortex breakdown 
may in fact occur in any vortex if the critical swirl and 
flow conditions are met. Hence the study of this phenomenon 
may have a practical application in any device containing a 
swirling fluid. 

Three distinct types of vortex breakdown have been re- 
ported by Sarpkaya [2]. These are the spiral, axisymmetric 
and double-helix types. Both the spiral and axisymmetric 
types initially form at a stagnation point along the vortex 
core, or centerline. These breakdowns mark the transition 
between laminar and “suse erate flows. 

In the case of spiral breakdown, the core appears to 
first suddenly stagnate, then kink to one side, and then 
form a spiral as the kink rotates. After a few revolutions 
the spiral becomes lost in the violently turbulent fluid 
flow. 

The axisymmetric breakdown also initially appears as 


a stagnation in the flow, followed by a roughly right-angled 
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kink. However, rather than forming a spiral, the rotating 
kink forms the forward boundary of an egg-shaped body, which 
fills from the rear, exhibits significant internal backflow 
and turbulence, and then discharges through the downstream 
side of the body. Depending on the magnitudes of both the 
circulation and flow rate (i.e., Reynolds number) in the 
vortex, the axisymmetric body may be closed from as little 
as 50% of its projected length, = up to 95%. When the body 
is only partially closed, turbulence commences at the exit 
om the: bodynnandiimkansMoutit oy ithesoutier reaches of the siti 
laminar vortex in a very rough conical shape. When the 
axisymmetric "bubble" is almost fully closed, the exit core 
appears to re-form for one or two body lengths, and then 
either kink again and form a spiral as described above or 
proceed immediately to full turbulence in a smooth conical 
shape. 

The double-helix type of breakdown does not include a 
definite seataaencn point in the core; rather, a mild dis- 
turbance causes a slight slowing and swelling of the core. 
The core then shears into a tape-like helix formation. Its 
appearance is that of two ribbons emanating from a slight 


expansion of the core streakline, then entwining about each 


z This type is sometimes referred to as a "near- 
axismmetric" type [2]. 
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other. In the case of the double~helix breakdown, which 
generally occur at low Reynolds numbers (Re = 700 to 1500), 
there is no immediate transition from laminar to turbulent 
flow. It may therefore be considered by some to be more of 

a flow instability, significantly different in nature from 
the more common and relatively abrupt spiral and axisymmetric 
vortex breakdown forms [2]. 

Shortly after Peckham and Atkinson reported their dis- 
covery, numerous researchers attempted to develop a suitable 
model to describe and predict the occurence of vortex break- 
Gown. The experimental contributions of the various workers 
have been well documented by Robertson [3], and will not be 
pursued herein. Alongside the experimental efforts, several 
theoretical analyses were developed, but the divergence of 
approach and initial assumptions have created significantly 
different and conflicting theories to explain the vortex 
breakdown phenomenon. Recently Hall [4] presented an excel- 
lent summary of these studies and compared the various anal- 
yses, highlighting their major strengths and weaknesses. 

He grouped the various proposals into three categories, in 
which the basic ideas were as follows: 

1. The phenomenon is in some sense like the separation 

of a two-dimensional boundary layer (Gartshore 
1962, Hall 1967, and others). 


2. The phenomenon is a consequence of hydrodynamic 
instability (Ludwieg 1962, 1965, and others). 


3. The phenomenon depends in an essential way on 


the existance of a critical state (Squire 1960, 
Benjamin 1962; 1967, Bossel 1967, 1969). The 
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proposals in this category have equivalent 

definitions for the critical state but a 

Lance (Om mMnicerpRmevayLoOnsMOlmalusy Sa end falcance,. 
However, after detailed analysis of the theories, Hall con- 
cluded that none of the proposals was adequate to completely 
explain what a vortex breakdown was, and when or where it 
would occur. 

Interestingly enough, Sarpkaya [5] noted that while the 
axisymmetric vortex breakdown basically corresponded to the 
finite-transition critical-state theory of Benjamin [6], 
the double-helix and spiral types were more in accord with 
Ludwieg's instability theory [7]. Carrying this one step 
further, Mager [8,9] presented a mathemat lea analysis based 
on the works of Hall and Sarpkaya, which essentially com- 
bined the finite-transition and hydrodynamic instability 
theories. In essence he showed that each theory was basi- 
cally correct as far as it went, but that each was only a 
portion of a larger theory. He proposed that a complete 
analysis should not stop when the breakdown occurred, but 
rather carry through from the laminar vortex to the onset of 
turbulent mixing [9]. This proposal will be discussed further. 

From the foregoing, it should be obvious that much is 
still not known about the vortex breakdown phenomenon. Not 
only is there controversy over how or why the breakdown 
occurs, but also on the degree of importance of the various 
parameters; e.g., velocity profile, magnitude of swirl, Rey- 
nolds number and the effect of axial pressure gradients. 


Sarpkaya [2] has shown that the first three of these are 
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quite significant, at least in the study of breakdowns in 
water. Several researchers including Hall [4], Jones [10], 
Lambourne and Bryer [11], Mager [8] and Sarpkaya [2], have 
concluded that an adverse (positive) axial pressure gradient 
in the flow has an effect on the occurrence of breakdowns. 
Lowson [12], as a result of his analysis and experiments, 
reported that "the final position of the majority of the 
practical breakdowns is due primarily to the reaction of 
pressure gradient effects on an already established break- 
down." However, as yet, no one has definitely verified either 
theoretically or experimentally the role that the axial 
pressure gradient is presumed to play. 

The purpose of the present investigation then was two- 
fold.) Lhe first objective was to examine what, af any, 
effect adverse pressure gradient has on vortex breakdowns. 
This was) done) by using oydindrical test sections of varying 
degrees of divergence to provide different pressure gradients 
BOGE WVOrvex EhOw. Tle BSc Cond: ObyieC uae West Om Dice) talay, 
examine Mager's theoretical analysis of vortex breakdown, 
and to determine whether it offered any more promise as a 
vortex breakdown predictor than its predecessors. Considerable 
data were taken on vortex breakdown position and type over 
a wide range of flow DeRANE SERS. GIG this was used to both 
graphically demonstrate the effect of the adverse pressure 
gradient and provide data for a computer analysis of Mager's 
hypothesis as suitably adapted to the present experimental 


system. 
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Ii. DESCRIPTION OF THE EXPERIMENTAL APPARATUS 


The experimental apparatus consisted of a test tank, 
water supply reservoir, centrifugal pump, overhead reservoir, 
Aer box, two dye injection lines, vane-angle adjustment 
mechanism, thermometer, vortex-—breakdown position/swirl 
angle measuring device, two rotameters and necessary piping 
anda valves. Connected as, Shown tm  Fireurer ir: 

A scale drawing of the test tank is presented in Figure 
Ag bw was constructed of 3/4" Plexiglas plates, held to- 
gether by screws and sealed with silicon sealant. The re- 
movable top was made in two parts for ease of handling during 
removal for changing the test sections, and was fitted to 
ElNe lOOGIY \aligia 2) SalalL@oin icvlloloeie (Yasleeus — ibmmSigael sealeioalioyertss alia 


the test tank were as follows: 


1. Sectioned baffle-plate - This was made. of 3/4" 
Plexiglas, and drilled with 39 rows of holes 


Sueppea mrom LO at the boutom to) 5/1 OL on 
top. This provided a nearly uniform velocity 
distribution in the test tank. Fine mesh 
screens were installed on each side of the 
baffle-plate sections to reduce the turbulence 
of the entering fluid. 


2. Adjustable vane mechanism - An additional 3/4" 
Plexiglas plate was attached to the upstream 
Well OtGNCuves trucker me A mMOnabare  CaurecuUllare Teabmieys 
9" in diameter, was set into this plate. 
Thirty-two fixed pins were attached to the ring 
ajay 2) SHAMS WIPO Galio@vllere eumeehy, IIL 25) Cleeicoos 
apart. Each of these pins fitted into a slot 
in a streamlined blade, or vane. The vanes 
were pinned to the inner plate, and pivoted 
together about their fixed pins as the movable 
ring was rotated by a hand-screw mechanism. 

The indicator on the hand-screw mechanism was 
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calibrated to read the vane angle directly to 
PhemMeacesinehclsinaGepicce Tan! wasp Capiapakes Oj 
changing vane angle from plus to minus sixty 
degrees, thus varying the amount of swirl 
imparted to the vortex. Blade dimensions and 
geometry are shown in Figure 3. 


3. Cone and bell-mouth - A six inch diameter cone 
was attached to the additional plate described 
above, inside the circle of adjustable blades. 
This cone and its matching bell-mouth formed 
the entrance to the test section. The center- 
line was bored to provide an entrance for the core 
dye stream. An outer dye stream (to show the swirl 
angle) was introduced through an adjustable probe 
located in the flow stream between the core and 
bell-mouth. To ensure that the bell-mouth (and 
the connecting test section) were centered on 
the cone, simple hand-screw adjustments were 
provided for both horizontal and vertical 
centering. 


h, Divergent test section and outlet tube —- Four 


pairs of these were made in matched sets, machined 
to close tolerances to ensure that there would 

be no obstruction to smooth flow through the bell- 
mouth, test section and outlet tube. To test 

the effect of pressure gradient, tubes of varying 
angles of divergence were manufactured and in- 
stalled. (The inlet diameter om all test sections 
was identical to that of the bell-mouth outlet, 
WOW MleStusee ton wl i(LS=1)) sada tenet 

of 10 3/8", and an outlet diameter of 2 174", 
providing a (half) angle of divergence (y) of 
1.38°. Test section #2 was made by adding a five 
ineh extension to) IS-l) maintaining the same angile 
of divergence, thereby giving it an outlet di- 
ameter Ot ouel/iC ee Senin ashe met nn Oma lu anaiaG 

Vii sO moc wanGden TS un hiade allkemest mmo uls. sens 

and y = 3.25°. Each test section was installed 
with an outlet tube that matched the test sec- 
tion diameter. The outlet tube was mated with 
aemavchinemMozZ Leave lus harkend, to constrict 

the flew invouthe ay 172" outlet pipe teading 

to the rotameter bank. 


Exterior fittings on the test tank were the inlet box, 
thermometer to measure water temperature, and a 26" over- 


flow column to maintain constant tank pressure and absorb 
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flow fluctuations. A small device to measure breakdown 
position and swirl angle rolled along a section of model 
railroad track permanently attached to the top of the test 
tank. <A one inch deep Plexiglas sighting bar was suspen- 
ded from a wheeled, counter-weighted cart. The sighting 
bar also provided a suitable location to mount a fixed pro- 
tractor and rotatable straightedge to measure swirl angle. 

TOMKEeDUMe SVS eMac mulee Onl GlovuURbDaMces as Oss Hsoiker. 
air vents were provided at several locations, and the entire 
test tank was mounted on foam rubber. 

Flow of water through the system can readily be traced 
by referring to Figure 1. City water entered the lower 
reservoir, was pumped to the overhead reservoir, flowed by 
BravesnyaelmuOmune) elMmMlcit mbox. manGdwenweredn une sesic suiaMmle.. . Ike 
then passed through the bavPLecolates and screens, slowly 
traveled the length of the test tank, passed through the 
adjustable vanes, and between the cone and bell-mouth where 
Lb ROrned saelaminare VOrve x une water sthen entered: the 
diverging test section and exited the tank through the outlet 
tube. After passing through the outlet nozzle, the water 
was carried to either (or both) of the rotameters, and 
discharged to a drain. Rotameter ratings at 100% flow were 
MORO aniGe iO calelonis) sere Manel ese vece ively) and provided 
the capability of achieving Reynolds numbers within a range 
OmeRe = 500m rom lon OO0r 

A photograph of the system is presented in Figure 4. 


Out of sight are the pump, lower reservoir, inlet box and 
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the overhead reservoir. The large cement column to the 
rear of the system supports the overhead reservoir, located 


approximately fifteen feet above the test tank. 
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III. EXPERIMENTAL PROCEDURE 


The initial phase of Bini s investigation consisted of 
determining the location of vortex breakdowns in tubes of 
Various angles of divergence for identical flow conditions. 
the variables which) could be controlled were circulation, 
Reynolds number, and degree of divergence of the test sec-— 
tion. Circulation (2), a measure of the vortex strength, 
was taken as the product of the tangential component of the 
velocity of the fluid leaving the vanes and the circumfer- 
ence of a circle drawn through the vane tips. From conti- 


nuity and geometry, this resulted in the following equation 


L245 


T = an(R. -m cos $) V 


total °2" v 


where » is the vane angle (read directly from the indicator 
in degrees) and the other terms are calculated or defined 
as shown in Figure 3. Circulation was normalized by defining 


a circulation number as 


Q = T/U,D, 


which may be expressed as 


} 
TTD (Ry - m cos $¢) sin w 


mm Oran ny, 


e 
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A simple computer program was written to calculate the 
corresponding circulation number for each vane angle 
serous. 

Reynolds number (Re) was used to normalize the flow 


rate through the test section, 


UD, 


Vv 


Re = 





where D,. was the test section diameter at the start of 
divergence (1.755"), U, was uniform velocity, and v the 
kinematic viscosity of the working fluid, water. UL was 
determined from the flow rate measured with a calibrated 
rotameter. As the experiment was run using city water, v 
did not remain constant, but varied as a function of water 
temperature. This variation has been taken into considera- 
tion in the calculation of Reynolds number by using the 
appropriate value of kinematic viscosity for the temperature 
measured. 

The experimental effort was conducted in four parts 
(one for each test section), and each part was conducted in 
an identical manner. The test tank was filled very slowly 
PrOmmEvne overhead mreservOur im Order tO; minimrze) Turbulence 
in the test tank. Air vents, located at strategic points, 
were opened to eliminate any air pockets that might other- 
wise have caused disturbances in the breakdown position. 


After the test tank was filled, it was allowed to remain 
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idle for approximately one hour to permit the remaining 
small amount of turbulence to dissipate.- When the fluid 
was presumed still, the water supply was resumed and water 
run at a low flow rate through the entire system. 

Data runs were made by selecting a suitable vane angle 
and holding it constant while varying the flow rate through 
the test section by carefully setting the rotameter valve 
or valves. Vane angle setting varied from every two or 
three degrees for lower circulation numbers to five degrees 
for higher circulations. The total range of circulation 
Nunbersieror wal west ysect ons! warded fom.) = 10550. to (5). 00 
(> = 13° to 55°), although only a portion of this range 
eould be) used for each test section. 

For a fixed vane angle, the rotameter setting was in- 
creased in small increments, and conditions in the test 
section recorded. When a vortex breakdown was present, the 


following information was recorded: 


1. Vane angle (° - to nearest 0.2°) 

2. Flow rate (% - to nearest 0.5%) 

3. Water temperature (°F - to nearest 0.2°) 

NS Sigalseal angle? (Ces tommacance si ace) 

5. Breakdown type, or types if in hysteresis region 

6. Breakdown position «(inches from start of test 
section divergence - to nearest 0.1") 

2 


It would have been preferable to allow the water to 
remain in the test tank from day to day, but this would have 
permitted dissolved air to settle throughout the system and 
disturb the flow. 


3 Swirl angle was taken as the angle of the tangent to 
the outerrdye stream at the: nearest crossing of the core 
prior to the breakdown. It is considered by some to be an 
indication of a critical flow condition in the vortex [2], 
but it was not used in this investigation. 
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7. Breakdown diameter (inches - to nearest 1/16") 
8. Percent that breakdown was closed, if axisym- 
metric type (% - to nearest 5%) 
OP GeO moO tahoe: vy nuntOGssame reer. if axisymmetric 
type (% - to nearest 25%) 
10. Comments - anything peculiar to the current 
data point 

Breakdown and swirl visualization was accomplished by 
injecting dye into the vortex as previously described. 
Commercial food coloring at full strength was used ag HRS 
dye, and although it was known to be slightly denser than 
water, the effect of the density difference was not 
noticeable. 

The most time consuming portion of the experiment was 
marking the location of the vortex breakdown. As the break- 
down is a dynamic phenomenon, its position oscillated within 
a narrow range, even though circulation and flow rate re- 
mained constant. Ali data points taken exhibited a slight 
random fluctuation in the axial direction. The average 
distance traveled during these fluctuations was approximately 
HO hws wNSmexpertence was Gained. siviwasi found that by, 
following the slow motion of the breakdown with the position 
indicator for about three minutes and mentally noting the 
limits, there would be roughly a 95% probability that the 
breakdown would remain in a 0.4" range. (This was true for 
both large and small flow rates.) The midpoint of this 
range was taken to be the (average) breakdown position. 

This procedure was further complicated in the hysteresis 


region when two (or more) breakdown types were alternately 


present for the identical flow and circulation conditions. 
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When this occurred, it was necessary not only to determine 
the average position for each, but also the relative propor- 
tion of time that each was present. This latter step was 
necessary to facilitate the preparation of graphs in which 
average breakdown position for a given flow and circulation 
setting would be shown. 

Finally, after two complete sets of data runs were made, 
photographs were taken of various breakdowns over a wide 
range of flow conditions. These were taken on 2 1/4" x 
2 1/4" negatives in both color and black and white, using 
only the lights installed to the rear of, and under the 
test tank. Several of these photographs have been reproduced 
as Figures 5 through 13, and 19 through 23. 

Upon completion of data runs and photography, the test 
tank was opened and cleaned, and a different test section 
fitted and installed. As the test sections were partially 
completed in advance, this step was accomplished within a 


day or two. 
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IV. ANALYTICAL CONSIDERATIONS 


An exploratory approximate analysis of the effect of 
adverse pressure gradient on the location of vortex break- 
down consisted of the consideration of the equations of 
conservation of mass and momentum, the representation of 
the approximate velocity profiles by suitable polynomials, 
and the calculation of the distance at which the axial 
velocity gradient increases rapidly leading to a stagnation 
point. 

The equations of motion for axially symmetric, incom- 
pressible flow, with gradients in the axial direction much 
smaller than those in the radial direction, may be written 


as [8] 


(ore) me) 0 (1) 
ye 

ee te (2) 
(uvr®) + (wwe®) = etr*(v, - DI, 3) 
(uwr),, + (wr), = -GIP, + Fest) py, (4) 


where Re. = Q R/V and all velocities are normalized by 
the maximum velocity On all pressures normalized by the 


maximum dynamic pressure 0.500. and all lengths normalized 
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by the radius of the tube at the start of the diverging 
test section R,- 
The above equations can be reduced to the following 


integral equations [8] 
Te =) 2/7Re (5) 


0 a ee 6 (SNe (alya) & © 


Tel, 22 1 (6) 
where 
6 
I, = J GW So wae Cho, 
re) 
ti 8 Pon & we a 
ma 7 CR = te eee, 
: (7) 
tis € 7 yea (=) the dr, 
8 9 
I55 Be eevee Cite. 
fe) 
and 
20s Gm 1 owe) 2 EQ)eat a dips tr) 
Oo fo) a Reo msIaaO© 
(8) 
with 
62 
APTS De mPa 20 dh a (9) 


where a = 52 and S = Wa ie Q@ are the normalized core area 


and the normalized circulation respectively. 
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The foregoing equations may be integrated further only 
if one assumes a set of suitable velocity profiles. This 
procedure is quite similar to that commonly used in employing 
the integral-momentum equation of senna One may do so 


by setting 


Weller eh Ob aieyy 38 A V(f, + Bf2), 


(10) 
v= 1,/6 


where 


n= r/6, eau = one a) 
(11) 


a) 
if] 


fe 2 = n°) 30, n)® 

These velocity profiles satisfy the boundary conditions at 

r= 0 and r = 6, and reduce to the swirling-potential-flow 
velocity distribution outside the vortex core. The parameters 
a and B are initially unknown and will be determined later. 
Finally, to. complete the formulation of the analysis, one 
needs to express the conservation of massflow within the 


diverging test section by writing 
MES SWAT = 22m) 4 Gp) 


“leva an example of this procedure, see Schlicting [13]. 
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where A is given by 
2 
e=) Ge aterniuny) aa (13) 


GUS tnevackalmGistance, wand yersevhe hali—angelke sor test 
section divergence in Equation (13). Evaluating the inte- 
grals given by Equation (7) by the use of the velocity pro- 
files, Equations (10) and (11), and inserting the various 
fy values in Equations (5,6,8,9, and 12), produces the 
following four equations: 

[w(K 4 PGK oF Ky 50°) - T?/Yala! +[aw°(K,, + QGHK 


ll Villon! 


12 
+ Pure (b db PRD.) IBY sb DEW + K..qg + K a) 
% ale Nae eae Oe Me temtalii ey vee 


+ 0.1 aW(1 - a) J]W' = [0] (14) 


[-0.2W(1 - a)Ja' + [0.2aWJa' + [A - 0.2a(1 - a) ]W' 


= [-WA'] (15) 


[(t/a)?a + b + d518 + By Re + Paci Te 
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+ [-(b5, + 28055) T°/ale' + [20°w - awlw' 


= [48(1 - a)/(a Re)] (16) 


[W(K, + aK,, ) Ja! + [aWK,.]la' + [aW(b 25 + ab), ) 18" 


21 
+ [a(Ky, + oK,,)] = [2/Re] (17) 
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The coefficients K and b 


fa 1,3 are contained in Appendix 


A. 
The above set of equations were written into a computer 
DBOrTam  andusolved™ samulitancousmy mien raat. ob! svande Ww". 


and then for a, a, 8, and W by writing, for example, 


with the initial conditions 


This process was continued until a singular point was reached, 
forcing the coefficient matrix of Equations (14) through (17) 
to zero. This condition occurs when a stagnation point is 
approached along the axis of the test section, and corresponds 
to that of vortex breakdown position [8,9]. 

It should be noted prior to a discussion of the utiliza- 
tion of the computer program that Mager attempted to calcu- 
daveline effect of swirl on nozzle flows through the use 
of a set of velocity profiles other than those used herein. 

An initial effort was made to use his formulas for the pre- 
diction of the vortex breakdown aloe aeaonk but it was discov- 
ered that there were not only several errors in his final 


5 


formulation, but the velocity profiles chosen by him did 


y Two typographical errors were found: (1) on page 
651, Table 1, coefficient C3 i should be 2Wa/5 vice 2WA/5. 
(2) On page 654, Appendix, function f should include the 


number 2697, vice 26977. a 
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not lead to a point of singularity even after the errors 
were corrected. Thus the analysis presented above, although 
similar to that pioneered by Mager, differs in many respects 
and certainly from the objective and results discussed 
by Mager in reference [8]. 

The computer program written from the above analysis 
has been included in a separate section following Appendix 
B. It was written in FORTRAN (IV), and can be most efficiently 
used on a time-sharing system communications terminal. 
Utilization of the program proceeded in the following man- 
ner. First the angle of divergence for a particular test 
section was entered into the program. Then the program was 
compiled and executed. The program was written so that it 
WOULGmDaAUSconGume quest n tMchmidaiveralls Wallies: OtuiCH brs) Olan ne). 
and eT: and the print spacing format; the latter term included 
so that the procedure could be watched iteration by iteration, 
or by using a large enough number, would print only values 
at the location of core stagnation. and Re were entered 
as the parameters of the particular breakdown position in 
question. The initial values of a and 8 were unknowns. 
They were found for this system by matching predicted break- 
down locations with those obtained experimentally for the 
same parameters. The procedure was to select a point, insert 
Phem conresponGinigno. wR anG minum On aac lO) ard. 
increase 8B until the maximum value of Z was obtained. Then 
using that 8, plus the previous 6, Re, and 2, &® was increased 


or decreased until the predicted value of Z corresponded to 
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that found experimentally. These computed values of a and 

B were used to predict breakdown locations for various values 
of 6, Re, and &. As will be discussed later, the foregoing 
exploratory analysis is capable of accurately predicting the 
effect of the increasing adverse pressure gradient and 
Reynolds number. This was all that was expected from the 
analysis since no precise calculation could be carried out 
without knowing the initial velocity profiles. The measure- 


ment of the latter is a prohibitively complex problem. 
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V. PRESENTATION AND DISCUSSION OF RESULTS 


A. EXPERIMENTAL RESULTS 

Photographs of the three primary types of vortex break- 
Gowns, axisymmetric, spiral and double-helix, are presented 
in Figures 5 through 7. These breakdowns could be produced 
in each of the four test sections, although the occurrence 
of the double-helix type was much less common in the test 
sections with higher angles of divergence, TS-3 and TS-4. 

In addition to the three major types, several sub-types 
of breakdowns and core disturbances were also noted. The 
first of these to be discussed appeared as an instantaneous 
Siaeniateon Om mbhe) Come low. tits: OCeunredvas a) low ke lopl 
phenomenon, and was quite unsteady. It simply appeared as 
a temporary slowing of the core that was immediately washed 
out before a breakdown could occur. As Reynolds number was 
increased with circulation held constant, ne next type that 
occurred was an "open bubble," shown in Figure 8. It formed 
from a stagnation point and exhibited the initial kink found 
in both the spiral and axisymmetric types. It differed, 
however, from both of these primary types as it exhibited 
some backflow behind the kink,. but did not possess the energy 
required to fully close its forward wall as in the axisymmetric 
case. Due to insuriiceient circulation, the core dye washed 
out of the breakdown after one-half revolution of the ro- 


tating kink. A further increase in Reynolds number produced 
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the spiral and near-axisymmetric types in the hysteresis 
region. These phenomena are demonstrated in Figures 9(a) 
Brae! Co) 6 Figure 9(a) is a spiral type and 9(b) a near- 
axisymmetric type. These two photographs were taken at the 
same flow and vane angle settings; the photographs were 
taken about one minute apart. It should be noted that in 
the lower range of Reynolds numbers the spiral form pre- 
dominated, while the axisymmetric form was the more common 
and persistent of the pair for higher flow rates. 

The hysteresis effect appeared to be caused by the 
streamlined shape of the vortex breakdown, and the action 
in the test section during hysteresis can be shown by the 
following cyclical pattern. Assuming that the breakdown is 
first a spiral type, it will move slowly forward (toward the 
start of tube divergence) and backwards in the axial direc- 
tion. When it moves forward beyond the limits of its average 
position (about + 0.2"), backflow will suddenly start and 
the axisymmetric type will form. The now streamlined shape 
will permit the pmeaidown to move forward as the fluid 
accelerates around it. It will steady in a new average 
position approximately 1" to 2" forward of that for the 
corresponding spiral type. As the near-axisymmetric type 
slowly expands and contracts, it will respectively move 
Slowly forward and backward in the axial direction. For 
some unexplainable reason, it will occasionally move back 
too far (out of its average range), and the breakdown will 


dissipate. When this occurs, the kink will continue to 
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rotate at the same speed, but the breakdown position in 

the test section will now permit only the formation of the 
spiral type. Then after a few moments of axial fluctuation, 
the spiral will move forward, re-form as an axisymmetric 
type, and the cycle will repeat. A similar hysteresis action 
was found to occur at very high Reynolds numbers for the 
near-axisymmetric and closed axisymmetric types. 

Still further increases in Reynolds number (to 
Re z 7000) produced the tightly closed axisymmetric type, 
shown previously as Figure 5. It has been described quite 
comprehensively in Ref. 2, and thus a description of this 
phenomenon will not be persued further herein. 

The above breakdown type transition pattern from 
instantaneous stagnation to the almost completely closed 
axisymmetric type was found to occur in all test sections. 
For the less divergent test sections (TS-1 and TS-2) at 
large circulation numbers, a completely different breakdown- 
transition pattern occurred. At extremely low Reynolds 
numbers (re z~ 800 to 1200) a core instability analogous to 
sine Pa eTaeOUe stagnation occurred, but it did not involve 
a stagnation point. As shown in Figure 10, the core streak- 
line was deflected from the test section centerline, sheared 
slightly, then re-formed and once again entered the test 
section centerline. This offset portion of the core did 
not rotate, but rather remained stationary in the position 
shown. This disturbance was accompanied by a large increase 


in swirl angle, shown in the photograph. Notice that 
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turbulence did not occur in the test section downstream of 
the disturbance. It had originally been assumed that this 
disturbance was simply a double-helix type for which the 

dye preferred to illuminate only one tail; however, even 

when great quantities of dye were admitted to the core, no 
evidence of an obscured double-helix type was found. A 
double-helix type could be formed from the core deflection 

by increasing the Reynolds number. When this was done, the 
slight swelling of the core occurred roughly at the point 

of the initial deflection described above. As in the core 
deflection case, the double-helix did not rotate - it appeared 
completely stationary. Although fluid was moving slowly 
through the test section, there was no visible fluid motion. 
A further increase in Reynolds number at constant circulation 
caused the double-helix type to re-form into a rather flat 
near-axisymmetric type (Figure 11). This breakdown also 
appeared motionless, except for very slow fluid flow within 
the "bubble" area behind the stagnation poinitn, anda sient 
wavering of the two tails emanating from the breakdown. 

OAes Casali sete breakdown did not rotate. A further increase 
in Reynolds number either caused a transition directly to a 
spiral or near axisymmetric breakdown, or in some instances 
at rather high circulation numbers (2 = 1.65 to 2.65), the 
formation of a rotating two-tailed axisymmetric type. 

Figures 12(a) and 12(b) show this type. The two photographs 
were taken seconds apart, and are shown to demonstrate the 


rotation of the tails. In this instance, the tails rotated 
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at approximately 39 rpm. Further increases in Reynolds number 
caused a transition to the near-axisymmetric type as noted 
above. Only very rarely was there sufficient circulation 

and flow strength to produce a completely closed axisymmetric 
bubble for high circulation breakdowns in test sections with 
low angles of divergence. This is probably due to the fact 
that for high circulation and Reynolds number combinations, 

an axisymmetric breakdown occurs extremely close to the bell- 
mouth eons (forward of the start of test-section divergence), 
where the vortex is not completely formed. 

One other breakdown form, shown in Figure 13, occurred 
for very high circulation numbers. It appeared to be a rather 
large axisymmetric type, but rather than filling from its 
closed downstream end, it filled from an axial, circumferencial 
ring located approximately two-thirds of the way back from 
the stagnation point. The nearly stagnant fluid inside the 
body completed the closed downstream portion of the breakdown. 
It emptied from the interior of the axial ring, rather than 
from the rearmost point of the body as would normally have 
been expected. No explanation for this phenomenon is 
Oishereds 

Figures 14 through 17 show curves of the average 
normalized breakdown position as a function of Reynolds 
number stor constant, circulation, -for.ach of the four test 
sections. Each curve represents the average of two data 
runs, weighted to account for the preference for one of the 


two (or three) breakdown types in the hysteresis region. 
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In addition, random curves have been selected, and for these 
all data points from which the average position curves were 
drawn are shown. These points demonstrate the scatter of 
the data for the two runs, the range and amount of hysteresis 
for different breakdown types, and finally the transition 
through the breakdown sub-types as Reynolds number was 
changed. 
Figure 18 graphically demonstrates the significant 
effect of adverse pressure gradient on vortex breakdown. 
In this figure the curves show average normalized breakdown 
position for various Reynolds numbers for each of the four 
test sections for a single circulation number, Q = 0.67. 
The significance of the curves for TS-2 and TS-3 is obvious - 
as adverse pressure gradient increases, the breakdown occurs 
closer to the origin of divergence for the same Reynolds and 
circulation numbers. However, the meaning of the locations 
of the curves for TS-1 and TS-4 are not so obvious. From 
the above, it would have been expected that the curves for 
TS-1 and TS-2 would coincide as they have the same angle of 
Hiyetaenes. and that the curve for TS-4 would be far to the 
left of that for TS-3 since the angle of divergence is so 
much greater. These expectations obviously did not occur. 
Inasmuch as TS-1l was five inches shorter than TS-2, 
the transition from the adverse pressure gradient caused 
by tube divergence to the nearly zero pressure gradient 
of the following straight section had the effect of slightly 
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adverse pressure gradient acting on the breakdowns in TS-1 
was ponennee Jess than that acting on those’ of Ts—2. 

The location of the curves for TS-4 relative to those 
Otel owas aveenestaquate!  Daiaulimp. = li awals) Oily, aint ex 
the outer dye injector was used to illuminate the flow 
pattern in the test section boundary layer that the actual 
fluid action in the tube was discovered. Figure 19 shows 
that boundary layer separation was occuring in TS-4. Rather 
than the boundary layer being confined to a very thin layer 
of stagnant fluid along the tube wall, there was a much 
thicker layer of fluid between the tube wall and the vortex. 
inches) thicker boundary, layer, reversed flow occurred, and 
this was actually the boundary layer that the vortex felt. 
The thin line that is visible in the upper portion of the 
tube in Figures 20 and 21 is the actual vortex boundary, and 
its angle of divergence is much less than that of the test 
section. Figures 22 and 23 show the fluid action in the 
boundary layer - Figure 23 being the more notable as it 
shows the early stages of the formation of a Taylor-Goertler 
vortex within that boundary layer. 

Figure 24 shows graphically the effect of the adverse 
pressure gradient on vortex breakdown. Here average norma- 
lized positions of breakdowns are plotted as a function of 
the test section divergence for test sections 2, 3, and 4, 
for various Reynolds and circulation number combinations. 


Notice that the breakdown positions for TS-1l have been 


36 

















ti 
@& 





ear At ienwotniootd add, do. ganton 3 


S~er Xo saontd’ eo gettos seri? naan 
geont of¢. ovissiet Rage “cl aevws edt 20, ose 







s9dte vino sew at ankitisd ei2up Sati? pr 


or 


worl? adt. gtsafiw lit ot beau anw totoetat ub 


























Teotés: od? Yee repel” eaabralen noliose drat edt ot x 
swore Of emg .berevooerh Baw “du oft at nota 
ender hab? nt ghinwose caw noltatsqee seysl yrebaee 
— teyst ntds Yar s od Benktnos prked novel caine 
dood 5 saw ated? flaw sdud- add _nole bivold oe 

a -AS2g0v of% boa Liew odud aft neowsed Biytt to age 1 
bas ,betz550 wolt beoetéver ,tevel yrabaued teNbs 

i <3fo% xotuot saa Jad? «ogat vrebruved ond eliausia 
7 . aad % mnbiqot oan m¢ at eldtatv eb gadtd ents mbt 
bas .winbevod réddev Intoe edt er [9 bnew OS etagee ‘ 


cu 


ined eto To fen? matt sweet doum af obiegtevid 36/88 


J | 
PS ei? af meléce bint? etd wore -€$ One-S? wag ae 
3? aé& oldsgvorm e#e1c ote gnted FS emglt = sewed 
sattiped-TolyaT 6 26 nottanrrol sit 20 asgede> eter ott | 
sayst uiaboiied ait? meta itw 
wha #02 -%6 wid yLisoltdqess awoke th SineaE 
“sic 4087 To!) wobMeotd xelmev po Jeeliers oe 
Gok ton i@ botdolg ots tnwobasstd Yo ened2taog beg 
Brot no) ot eohegtev th nolzaee 8982 
olsenidua® sed teclinortes bra ebfoagek egoitsyv 4607 


f he ny si stteoa mwdtbiseesse od: 





iti enttol 
a6 


included, and show that the effective angle of divergence 
for that test section was variable, but approximately equal 


to 1.30°, slightly less than its actual divergence, 1.38°. 


B. ANALYTICAL RESULTS 

The results of the exploratory analysis of vortex break- 
downs precscrLOnmslsmpreeSente duck here! 25.00 Uneiass) tel eure. 
several computer predicted breakdown locations are plotted 
with those obtained experimentally. Although there is sub- 
stantial deviation between predicted and experimental results, 
particularly at higher Reynolds numbers, the experimental 
curves are fairly well represented by the computer results. 
The primary reason for the deviations noted is attributed 
to the difficulty of exactly reproducing the initial velocity 
PrOmAMes soroueh chew use son muhem parameters sOelpiy) aia Of eara 
addition, the mathematical complexity of introducing the 
boundary layer development along the test section wall and 
the difficulty of assuming an appropriate eddy viscosity for 
the turbulent vortex core further reduce the accuracy of 
this analysis. Taken together, the approximations made 
herein preclude a more precise analysis based on the quasi- 
steady state momentum equations. Thus this analysis can 
only be considered an attempt to explore the general influence 
of the adverse pressure gradient, rather than a model to 
accurately predict vortex breakdown location. 

It is presumed that the precise measurement of the initial 


velocity profiles by more sophisticated techniques (possibly 
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through the use of laser-Doppler annemometry) would signifi- 
cantly improve this analysis, and thereby produce much more 


accurate results. 
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VI. CONCLUSIONS 


1. Adverse pressure gradient has a significant effect 
on vortex breakdown position. An increase in adverse pres-— 
sure gradient has the same effect on breakdown position as 
an increase in vortex circulation or mean flow rate —- that 
of shifting the breakdown location upstream, toward the 
origin of the vortex. 

2. An inerease in the angle of divergence of the test 
section beyond some specified point does not result in an 
increase in the impressed adverse pressure gradient. In- 
stead, the boundary layer separates and reversed flow occurs 
on the tube wall. This thicker boundary layer, vice the 
test section wall, then defines the outer extremity of the 
swirling flow, and thereby limits the effective adverse 
pressure gradient acting on the vortex. 

3. An increase in adverse pressure gradient does not 
directly affect the type of vortex breakdown created. It 
does, however, cause the constant circulation curves to 
shift down and to the left (toward the origin) on a plot 
of normalized breakdown position versus Reynolds number. 
Curves for higher circulation numbers are thus forced into 
the lower Reynolds number region of instability. This in 
turn restricts the formation of the relatively less stable 
breakdowns, particularly those afi the double-helix and 


related types. 
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4, In addition to the three primary vortex breakdown 
types (spiral, axisymmetric and double-helix), there exist 
numerous sub-types of breakdowns and Bore disturbances. 

5. The model developed to predict the location of 
vortex breakdowns performs only marginally well. The 
assumptions on which it is based, particularly that of the 
initial velocity profiles of the vortex, require further 
refinement. The method does, however, show promise, and 


provides a direction for further investigation of this 


matter. 
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Figure 4. Photograph of Experimental Apparatus 
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Figure 9(a). Example of a Spiral Breakdown in the Region 
of Hysteresis (TS-1, y = 1.38°, Re = 2100, 
Q = 1.04) 





Figure 9(b). Example of an Axisymmetric Breakdown in the 
Region of Hysteresis (TS-1, y = 1.38°, 
Re = 2100, 2 = 1.04) 
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Example of a "Flat Bubble" 


Figure 11. 


(TS-1, y = 1.38°, Re = 1500, 2 = 1.65) 
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Figure 12(a). Example of a "Two-Tailed Bubble" 
GUS 2k ry =) 538°. Rel= 1300, 9 = 2.66) 





Figure 12(b). Example of a "Two-Tailed Bubble" After 90° 
Cr howardon (US—2 5 ye= 38°.) Rey = 1,300) 
i=" 21616) 
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Figure 14. 


Breakdown Locations for Test Section #1] 
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Figure 15. Breakdown Locations for Test Section #2 
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Figure 17. Breakdown Locations for Test Section #4 
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Figure 18. Constant Circulation Number Curves of 
Breakdown Locations for Test Sections 
#1 Through #4 
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Figure 25. Actual Vortex Breakdown Position Versus 
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COMPLTER PROGRAM FOR PRECICTING VORTEX BREAKDOWN LOCATION 
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